*Table A.39: Meal and Snack treatment effects on “Added Sugar” - correction for attrition
cd ""

**NMESG
*Use the main dataset
use master_data_adults_child_s13456_shared.dta, clear
drop meal_control snack_control
gen meal_control=1 if treatment==2
replace meal_control=0 if treatment==1
gen snack_control=1 if treatment==3
replace snack_control=0 if treatment==1


drop if type_id=="X" // **Dropping newborns born during the three years of study
drop if id=="3885553B" //**error in the id number and empty data
drop if id=="3885553G" //**error in the id number and empty data
replace treatment=2 if id=="38855553B" //**retrieving the correct treatment for this id
replace treatment=2 if id=="38855553G" //**retrieving the correct treatment for this id
drop if id=="101aa22c" //**this id does not exist

*
ge food15_other_veg_new=food15_other_veg
replace food15_other_veg_new=0 if food15_other_veg==. & energykcal!=.

ge food16_fruit_new=food16_fruit
replace food16_fruit_new=0 if food16_fruit==. & energykcal!=.

su food15_other_veg 
su food16_fruit 
su food15_other_veg_new 
su food16_fruit_new

*
drop portion_veg_wisp portion_fruit_wisp
gen portion_veg_wisp=food15_other_veg/80
gen portion_fruit_wisp=food16_fruit/80
gen fuit_and_veg_g=total_veg_fruit_eatwell if session==1
replace fuit_and_veg_g=total_veg_fruit_eatwell if session==3
replace fuit_and_veg_g=food15_other_veg+food16_fruit if session==4

********************************************************************************
ge meal=(treatment==2)
replace meal=. if treatment==.

ge snack=(treatment==3)
replace snack=. if treatment==.

ge session1=(session==1)
replace session1=. if session==.

ge session3=(session==3)
replace session3=. if session==.

ge session4=(session==4)
replace session4=. if session==.

ge session5=(session==5)
replace session5=. if session==.

ge session6=(session==6)
replace session6=. if session==.

ge meal_session1=meal*session1
ge meal_session3=meal*session3
ge meal_session4=meal*session4
ge meal_session5=meal*session5
ge meal_session6=meal*session6


ge snack_session1=snack*session1
ge snack_session3=snack*session3
ge snack_session4=snack*session4
ge snack_session5=snack*session5
ge snack_session6=snack*session6


set more off
********************************************************************************
xtset id_nb session
***Panel A: Children
**Column 1 Energy 
*Results
ge FV=food15_other_veg_new+food16_fruit_new

xtreg nmesg session3 session4 session5 session6 meal_session3 meal_session4 meal_session5 meal_session6 snack_session3 snack_session4 snack_session5 snack_session6 if child==1, cluster(hh_id) fe

ge t=session
replace t=t-1 if session>2
ta t

keep if child==1
sort id_nb t
by id_nb: generate n2 = _N

xtset id_nb t
ge x=t==5
egen lastsession=max(x), by(id_nb)

keep id_nb treatment session session3 session4 session5 session6 meal_session3 meal_session4 meal_session5 meal_session6 snack_session3 snack_session4 snack_session5 snack_session6 child hh_id nmesg lastsession n2 t
save temp_before_trimming_nmesg, replace


use temp_before_trimming_nmesg, replace
merge m:1 id_nb using trimmed_samples_s3DD_AS
destring lower_trimDD upper_trimDD, replace

xtreg nmesg session3 session4 session5 session6 meal_session3 meal_session4 meal_session5 meal_session6 snack_session3 snack_session4 snack_session5 snack_session6 if child==1 & (session==1 |  session==3), cluster(hh_id) fe
xtreg nmesg session3 session4 session5 session6 meal_session3 meal_session4 meal_session5 meal_session6 snack_session3 snack_session4 snack_session5 snack_session6 if child==1 & lower_trimDD==1  & (session==1 | session==3), cluster(hh_id) fe
xtreg nmesg session3 session4 session5 session6 meal_session3 meal_session4 meal_session5 meal_session6 snack_session3 snack_session4 snack_session5 snack_session6 if child==1 & upper_trimDD==1 & (session==1 | session==3), cluster(hh_id) fe


use temp_before_trimming_nmesg_v2, replace
merge m:1 id_nb using trimmed_samples_s4DD_AS
destring lower_trimDD upper_trimDD, replace

xtreg nmesg session3 session4 session5 session6 meal_session3 meal_session4 meal_session5 meal_session6 snack_session3 snack_session4 snack_session5 snack_session6 if child==1 & (session==1 |  session==4), cluster(hh_id) fe
xtreg nmesg session3 session4 session5 session6 meal_session3 meal_session4 meal_session5 meal_session6 snack_session3 snack_session4 snack_session5 snack_session6 if child==1 & lower_trimDD==1  & (session==1 | session==4), cluster(hh_id) fe 
xtreg nmesg session3 session4 session5 session6 meal_session3 meal_session4 meal_session5 meal_session6 snack_session3 snack_session4 snack_session5 snack_session6 if child==1 & upper_trimDD==1 & (session==1 | session==4), cluster(hh_id) fe


use temp_before_trimming_nmesg_v2, replace
merge m:1 id_nb using trimmed_samples_s5DD_AS
destring lower_trimDD upper_trimDD, replace

xtreg nmesg session3 session4 session5 session6 meal_session3 meal_session4 meal_session5 meal_session6 snack_session3 snack_session4 snack_session5 snack_session6 if child==1 & (session==1 |  session==5), cluster(hh_id) fe
xtreg nmesg session3 session4 session5 session6 meal_session3 meal_session4 meal_session5 meal_session6 snack_session3 snack_session4 snack_session5 snack_session6 if child==1 & lower_trimDD==1  & (session==1 | session==5), cluster(hh_id) fe
xtreg nmesg session3 session4 session5 session6 meal_session3 meal_session4 meal_session5 meal_session6 snack_session3 snack_session4 snack_session5 snack_session6 if child==1 & upper_trimDD==1 & (session==1 | session==5), cluster(hh_id) fe


use temp_before_trimming_nmesg_v2, replace
merge m:1 id_nb using trimmed_samples_s6DD_AS
destring lower_trimDD upper_trimDD, replace

xtreg nmesg session3 session4 session5 session6 meal_session3 meal_session4 meal_session5 meal_session6 snack_session3 snack_session4 snack_session5 snack_session6 if child==1 & (session==1 |  session==6), cluster(hh_id) fe
xtreg nmesg session3 session4 session5 session6 meal_session3 meal_session4 meal_session5 meal_session6 snack_session3 snack_session4 snack_session5 snack_session6 if child==1 & lower_trimDD==1  & (session==1 | session==6), cluster(hh_id) fe
xtreg nmesg session3 session4 session5 session6 meal_session3 meal_session4 meal_session5 meal_session6 snack_session3 snack_session4 snack_session5 snack_session6 if child==1 & upper_trimDD==1 & (session==1 | session==6), cluster(hh_id) fe
